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Abstract 

Decoherence of a shallow donor electron spin in silicon caused by electron-lattice 
interaction is studied. We find that there are two time scales associated with the evo- 
lution of the electron spin density matrix: the fast but incomplete decay due to the 
interaction with non-resonant phonons followed by slow relaxation resulting from spin 
flips accompanied by resonant phonon emission. We estimate both time scales as well 
as the magnitude of the initial drop of coherence for P donor in Si and argue that the 
approach used is suitable for evaluation of phonon induced decoherence for a general 
class of localized spin states in semiconductors. 
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Rapid development of quantum information science has revived interest in study 
of electron and nuclear spin polarization, specifically in condensed matter physics. Po- 
tentially powerful applications of electronic and nuclear spins in quantum computing 
[1] and spintronics [2] have recently been proposed. Among these quantum informa- 
tion devices the ones utilizing spins in semiconductors seem to be the most promising 
candidates for practical implementation. Indeed the existing semiconductor technology 
may someday allow for scalability of the elementary quantum logic gates into complex 
integrated circuits needed for information processing. 

The problem of fabrication of the semiconductor quantum processing devices opens 
up a number of experimental and theoretical challenges. These include control problems, 
such as manipulations by interactions of individual spins with external field and with 
each other [3] (in the present case the 1/2-spins represent bits of information and will be 
conventionally referred as qubits) , the initialization and readout of the quantum register 
[4], as well as precise positioning of the qubits (within several A) in a semiconductor 
host - the issue that requires significant advances made in contemporary lithographic 
technologies. 

Another requirement for implementation of semiconductor devices for quantum in- 
formation processing is a long time memory of the state of a quantum register, i.e. small 
decoherence of the spin qubits due to the interaction with environment. Several theo- 
retical works that study decoherence in systems contemplated for quantum information 
processing have been made recently [5] . There is an outstanding set of experiments on 
spin relaxation processes in GaAs systems [2,6]. It was discovered that at certain doping 
concentration the transverse relaxation time of the electron spins injected into con- 
duction band by optical pumping techniques increases by several orders of magnitude 
and exceeds 100 ns under certain conditions. This remarkable phenomenon suggests 
possibilities of its use in the new generation of electronic nano-devices that utilize po- 
larized spins as information carriers and can become the building blocks of the future 
electronics. 

In the present work we focus on study of decoherence of localized electronic spins 
in semiconductors. The localized electronic or nuclear spins have been proposed to play 
role as qubits in several quantum computing architectures [7]. It has been long known 
[8-10] that in certain semiconductors (such as Si or Ge) the spins of electrons bound 
to donor atoms have very long longitudinal relaxation time T\ due to the inversion 
symmetry of the semiconductor host as well as due to the localized nature of their 
orbital states. This makes these electron spin states an excellent choice for qubits of a 
quantum computer. However, the transverse relaxation of the spin polarization in these 
systems, to the best knowledge of the authors of the work, has never been carefully 
investigated. 

We calculate the dynamics of the off-diagonal density matrix elements of the elec- 
tronic spin due to its interaction with acoustic phonons. These density matrix elements 
represent coherencies of the system intrinsic for a quantum mechanical superposition 
state, where individual states enter with certain weights/probabilities and phase shifts. 
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Interaction with environment, i.e. phonons, destroys the superposition by randomizing 
relative phases of the "up" and "down" states, thus leaving the state of a qubit in a 
statistical mixture not suitable for quantum computing. 

We find that there are several time scales in the process of such decoherence. The 
smaller time scale is defined by the size of the orbital localized wave function of the 
electron (see below). The dephasing of this sort is governed by the interaction with 
the off-resonant phonons, primarily with wavelength comparable to the effective Bohr's 
radius of the donor. Decoherence on this time scale - the drop of the off-diagonal 
elements of the spin density matrix is fast, but incomplete, thus leaving the system in a 
slightly mixed state. We term this type of decoherence as phase uncertainty of the state 
due to the interaction with environment. The later dynamics, i.e., on a much larger time 
scale, is determined by the conventional relaxation processes associated with spin-flips 
produced by the resonant phonons studied in [9,10]. 

We consider the spin-lattice interaction Hamiltonian for a single electron donor in 
a semiconductor. Such Hamiltonian can be generally written as: 

7~L S -l = Afl B [(J x i XX H x + OylyyKy + a z e zz H z ] 

+ Bn B [exy((7xHy + (TyH x ) + c.p.}+CnBV HA, 

where hb is Bohr's magneton, cr's are Pauli matrices, e a p is the effective strain tensor 
averaged over the ground state effective mass function of the donor, A is the effective 
dilatation and cp stands for "cyclic permutation" [9]. Constants A, B and C depend 
on the specific properties of the host semiconductor as well as of the donor atom. For 
simplicity we assume that the magnetic field is parallel to the [001] crystal axis, i.e., 
H x = H y = 0. In this case there is no spin relaxation due to the first term in Eq. 
(1). Such alignment of the magnetic field can be easily achieved experimentally. The 
anisotropy of the spin relaxation rate depending on the direction of the magnetic field 
for P donors in Si has been directly observed [8]. Then the Hamiltonian (1) is reduced 
to the following form (in interaction representation): 

Hs-i(t) = (JzQzit) + <r+Q-(t) + a-Q+(t) , (2) 

where a± = (a x ± ia y )/2 and 

Q z (t) = fi B H[Ae zz (t) + CA(t)} , (3) 

Q ± (t) = BfiBHe^'^it) T ie yz (t)} , (4) 

with u s = gfi B H, where g is the effective electronic g-factor of the host semiconductor. 
It is assumed here that the ^-factor that determines the Zeeman splitting for the donor 
electron is the same as that for a free electron. 

The effective strain tensor in (2)- (4) can be expanded in waves with definite wave 
vectors q and polarizations n p 

e a/3 = /(q,n p )e a/3 (q,n p ) , (5) 
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where e a/ g(q, n p ) can be written in terms of phonon creation/annihilation operators as 

e«/3(q,n p ) = ( — -5 ) (n pa qf3 + n p pq a )e~ lLJ ^ n p t al l:np + H.c. . (6) 

In Eq. (6), [a q>nj> , a q , >n ,J = <W<5n p ,nV w q,n p is the phonon frequency, p is the crystal 
mass density, V is the normalization volume, and H.c. denotes the Hermitian conjugate 
part. Function /(q) is the matrix element of the wave with definite q taken between the 
effective mass functions of the donor [10]. Only those lattice vibrations that have wave- 
length comparable or less than the spatial spread of the donor electron wave function 
contribute to the interaction (1). In the spherical donor ground state approximation 
/(q) = 1/[1 + a 2 B q 2 /4] 2 , where as is the effective Bohr radius for the donor [10]. Even 
though this is not the case for many semiconductors, such as Si or Ge, we will use the 
above form of /(q) for order of magnitude estimates. 

We are interested in the evolution of the reduced density matrix of the electronic 
spin, p = Tt p h (p s +ph), with trace taken over the phonon degrees of freedom. Here p s + P h 
is the density matrix of the full system whose evolution is governed by the interaction 
Hamiltonian (2)-(4). The temporal evolution of the reduced density matrix can be 
described by the Markovian master equations [11]: 

1 y 

P$) = Yl / dT [(<Ti<7jP(t) ~ <7jp(tWi)(Qi(t)Qj(T)) 

ij ° 

-{o l p{t)v J - pWajViXQjWQiit))] . (7) 

In (6) i,j = z, +, — , while ( ) indicates averages taken over phonon occupation numbers, 
wave vectors and polarizations. Equations (7) are generally valid on a timescale large 
compared to the correlation time of the phonon bath. However, it can be shown that 
when only the a z term is retained in the Hamiltonian (2), the decoherence, i.e. the decay 
of the off-diagonal elements of the spin density matrix, for the exactly solvable model 
[12] coincides with the solution obtained from (7). Thus we expect that the Markovian 
assumption built in (7) will still be applicable for the short time regime. Moreover, for 
the specific case of P donor in Si crystal considered below, the dimensionless constant 
A is an order of magnitude greater than B and therefore the mechanism of decoherence 
resulting from coupling to a z only can be considered separately. The master equation 
(7) allows us to combine both mechanisms, i.e. phase relaxation via the diagonal and 
off-diagonal spin components, into a single framework and so in what follows we will 
pursue with the solution of (7). 

In general, equation (7) represents four coupled linear differential equations. Their 
solution is a cumbersome analytical problem. These equations, however, are greatly 
simplified if one assumes that the semiconductor host medium is elastically isotropic. 
Such assumption is commonly used in evaluation of phonon related relaxation processes. 
For the isotropic phonons the only non vanishing correlators in (7) are (Q z (t)Q z (r)) , 
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(Q-(t)Q + (r)) and (Q + (t)Q-(r)). As a result, the equations for the off-diagonal ele- 
ments decouple from the equations for the diagonal density matrix elements. For the 
off-diagonal element P12 one easily obtains that 



with 



Pi 2 (*) = Pi2(0)exp[-r(£)] (8) 

t t' 

T(t) [ dtf f dr (2{[Q z (t'), Q z (r)] + ) + ([Q_(f ), Q+(r)]+» • (9) 

n Jo Jo 

We evaluate correlators entering (9) for P donors in Si. It was shown by Roth [9] that 
there are two major mechanisms contributing to the donor electron spin relaxation and 
thus determining the interaction Hamiltonian (1),(2): g-factor modulation due to inter- 
valley mixing caused by the phonon induced strain as well as (7-factor modulation in each 
valley [9]. The first mechanism determines constants A and C in (1),(2). The estimates 
of [9,10] for Si:P give A 2 ~ 10 and C ~ —A/3. The second mechanism is responsible 
for Hamiltonian with coefficient B in (1), (2). The estimate in [9] gives B 2 ~ 10 _1 . 
While the first term clearly dominates at short times, it saturates for longer times and 
so the second mechanism prevails. Therefore the first two correlators in (9) that are 
proportional to A 2 are evaluated exactly, while the last two terms (with coefficient B 2 ) 
are evaluated in secular approximation (with only spin flips contributing) that is valid 
only for long times [11]. After performing averaging over polarizations and directions 
of phonon wave vector q one obtains that the real part of T(t) is given by 

Rer(t)~r o7 ,(t) + — , (10) 



where 



4A 2 p 2 B H 2 ^ 
lh-K 2 a 2 B pc^Ti ' 

lM = a \ jT dqqf( q ) (1 - cos (c t tq)) coth 



(ii) 



. v , ; dqqf z (q) (1 - cos (citq)) coth ^ ) , (12) 



and 



AB 2 (ji B H)' 

07TM p 



r- 1 



1 2 



cothgg). ,13) 



In (10)-(13) c\ and Ct are longitudinal and transverse sound velocities, T is the temper- 
ature. For silicon c\ > Ct and so the last terms in (12) and (13) can be neglected for an 
order of magnitude estimate. The integrals (12) converge due the presence of /(g) which 
imposes a natural cut-off q < a^ 1 on the magnitude of the wave vector q. At the same 
time in evaluation of r^ 1 given by (13), f(q) does not contribute because only resonance 
phonons conditioned by Tic t q res = 2psH enter in (13), and f(q res ) — 1- Note that in 
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Eqs. (11) and (13) we also assumed that the average fir-factor is approximately 2, which 
is the case for Si. It also should be emphasized that at low temperature (for quantum 
computing applications we assume that T ~ 10 mK) coth (hciq/kT) ~ 1 over the range 
of integration < q < . Then it can be seen that the function 7 2 (t) approaches its 
maximum value equal 1/3 for t > a B /ct ~ 10 -12 s, meaning that the process is very 
fast and therefore in an applied magnetic field the phase uncertainty of order To/3 is 
practically always present. A similar fast, but incomplete decoherence was found in 
the study of phase relaxation of nuclear spins in quantum Hall ferromagnet [13]. We 
argue that this phase uncertainty generated by the first term in (11) is common for 
most phonon-assisted spin decoherence mechanisms and should be accounted for when 
evaluating whether a system is a good candidate for quantum computing. 

For H = 3 T, the magnetic field typical for quantum computing architectures [3,7], 
we find for P donor in Si that To/3 ~ 10~ 9 and tr ~ 1 s; see Fig. 1. Note that our Eq. 
(13) for the conventional spin relaxation has a different from Refs. [9,10] dependence 
on the applied magnetic field, H 5 instead of H in [9,10], which comes from the fact 
that we are interested in a very low, effectively zero temperature regime. Indeed, for 
high temperatures in (13), coth {p B H /k B T) ~ (k B T / p B H) and we recover results of 
[9,10]. Our results indicate that the P donors in Si are excellent candidates for qubits 
of a quantum computer as both the phase uncertainty and the spin-flip rate satisfy the 
10 -5 "clock speed to decoherence" criterion [14] for reliable quantum computing. 

Another aspect not considered so far is decoherence due to the interaction of the 
donor electron spin with optical phonons. Hamiltonian (1) derived in [9, 10] does not 
include optical phonons as they do not contribute to the real, i.e., spin flip involving 
relaxation processes, due to the presence of the gap uq in their spectrum (uq is optical 
phonon frequency at wave vector and ujq ^> HbH). Therefore coupling of the optical 
phonons to x and y components of the spin will create virtually no decoherence. In 
order to estimate decoherence produced by the optical phonons, we derive the effective 
coupling of the donor electron spin with optical deformation. The procedure is identical 
to the derivation of the Hamiltonian (1) for interaction with acoustic phonons, that 
was obtained perturbatively in [9,10]. Then we obtain that (see [10] for the details of 
identical derivation) coupling of a z to optical phonons is given in terms of the interaction 
Hamiltonian 

n opt = Dfi B Ha z ^^-, (14) 

where dimensionless constant D is of order (gi—gt), the difference between the longitudi- 
nal and transverse to the field g-factors, D ~ 5 x 10 -3 for Si, (W op t) is the matrix element 
of the optical deformation potential taken between donor electron effective mass ground 
state functions and A is the energy splitting between the orbital singlet ground state 
and the excited doublet states of the donor, A ~ 10~ 2 eV [15]. The decoherence due to 
the interaction (14) can be evaluated using Eqs. (9)-(12). One obtains that decoherence 
due to Hamiltonian (14) is also incomplete with r^* ~ [(gi — g t )p B H A] 2 /[p A 2 a b uj q?i 3 ], 
where A is the optical deformation potential constant. For typical parameters of Si one 
gets that Tq P * ~ 10 -11 and so the optical phonons create effectively no decoherence for 
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the localized electron spins in silicon. 

One should note that the case of P donor in Si is unique even among the semicon- 
ductor crystals with inversion symmetry of the lattice. Indeed the spin-orbit interaction, 
that mediates interaction of the spin with lattice vibrations, in Ge is several orders of 
magnitude stronger than in Si. As a consequence the effective coupling constants A and 
B in the Hamiltonians (1), (2) increase by several orders of magnitude [9,10], resulting 
in a much higher, presumably by a factor 10 6 , relaxation rate and initial decoherence. 
Finally, it should be pointed out that the nuclear spins can significantly contribute to 
decoherence of localized electronic spins by creating fluctuating magnetic fields via the 
diffusion of the nuclear spin polarization due to nuclear dipole-dipole coupling [16]. Such 
mechanism, though ineffective for longitudinal spin relaxation (Ti) because of large dif- 
ference in electron and nuclear Zeeman energies, can be important in phase relaxation 
processes, and should be subjected to detailed study. 
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Figure Caption 



FIG. 1. Schematic dependence of the spin off-diagonal density matrix element on 
time for P donor in Si. The initial, nonexponential coherence drop on a picosecond scale 
is followed by exponential, relaxation driven decay. 



